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The ratio of the I2 C to C isotope is a good tracer of the Galaxy (e g PaUa et al their Fig 4 and 5 } Such models 
amount of stellar processing in low and intermediate mass use the assumption that mixing in the ISM is comp iete and re- 
stars during the asymptotic giant branch (AGB) phase. The stricted tQ material at a giyen galactocen tric distance (azimuthal 
third dredge-up increases I2 C/ I3 C to about ~ 300, while stars mixing) A comparison of results between our Galaxy and other 
massive enough to undergo "hot bottom burning" bring their galaxies win giye important data for the nuclear process ing his- 
surface carbon isotopic ratio to CNO processing equilibrium, tory ^ e g ^ Tosi 200a Prantzos 20 01). 
12 C/ 13 C~3. The carbon yields of AGB stars are sensitive func- There are a large number of measur ements of the 12 C/ 13 C 
tions of their initial masses and metallicities (e.g. Renzim & ratio from radio astronomy data (see , e.g., Wilson & Rood 
VoliU981] or Casassus & Roche |2001] for a population synthe- rrggfl However, these results may be affected by interstel- 
sis approach and references to model results). lar chemistry in two ways . These are chemical fractionation 
A measurement of the I2 C/ I3 C ratio in the interstellar (which enriches molecules in 13 C and thus lowers the ratio) and 
medium (ISM) gives important data for the total amount of low selective dissociation (which destroys the rarer species more, 
mass stellar evolution and subsequent enrichment of the ISM an d thus raises the ratio). The 13 C enrichment is due to frac- 
in the Ga laxy T he terrestrial value of 12 C/ I3 C is 90 (Rosman donation of the CO molecule via the ion-molecule process, 
& Taylor 1998 1. The 12 C/ 13 C ratio is a cornerstone of models 

of the nuclear history of our galactic ISM, that is, the Galactic 13 C + + 12 CO ^ 12 C + + 13 CO + AE, (1) 

Chemical Evolution. Models of chemical evolution predict a • , , T , „ m -? i • , ,,• 

F with AE/k = 35 K and k = 2 10" lu cm 3 s" 1 in both directions 

(Watson et al. 1976 1. If the reaction reaches equilibrium, 
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Abstract. The 12 C/ 13 C isotope ratio is a tracer of stellar yields and the efficiency of mixing in the ISM. I2 CH + / 13 CH + is not 
affected by interstellar chemistry, and is the most secure way of measuring 12 C/ 13 C in the diffuse ISM. R =' 2 C/ 13 C is 90 
in the solar system. Previous measurements of 12 CH + /1/13957. 7,4232. 3 and I3 CH + /1/13958. 2,4232.0 absorption toward nearby 
stars indicate some variations in 12 C/ 13 C, with values ranging from 40 to 90 suggesting inefficient mixing. Except for the 
cloud toward f Oph, these R values are strongly affected by noise. With UVES on the VLT we have improved on the previous 
interstellar 12 C/ 13 C measurements. The weighted 12 C/ 13 C ratio in the local ISM is 78.27 ± 1.83, while the weighted dispersion 
of our measurements is 12.7, giving a 6.9 cr scatter. Thus we report on a 6.9 cr detection of 16.2% root-mean-square variations 
in the carbon isotopic ratio on scales of ~100 pc: R = 74.7 ± 2.3 in the £Oph cloud, while R = 88.6 ± 3.0 toward HD152235 in 
the Lupus clouds, R = 62.2 ± 5.3 towards HD1 10432 in the Coalsack, and R = 98.9 ± 10.1 toward HD170740. The observed 
variations in 13 C/ 12 C are the first significant detection of chemical heterogeneity in the local ISM. 

Key words. ISM: abundances - ISM: clouds — ISM: molecules ISM: Coalsack - ISM: Lupus cloud - ISM: <fOph cloud 

1 . Introduction decrease in the 12 C/ 13 C ratio with time for a given galactocen- 

tric distance, and a decrease with galactocentric distance in the 
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In cold gas, with T % 35 K, CO will be enriched in 13 C com- 
pared to C + . The effect is observed in CO absorption towards 
£Oph (Sheffer et al. [TffHl . pOph and ^fOph (Federman et al. 
2003 1. Thus molecules whose synthesis involve CO at a later 
stage than C + will show the CO fractionation, and vice-versa 
(e.g. Watson 1978 1. The corresponding 13 C depletion in C + is 
diluted in environments where n(C + ) » n(CO), such as diffuse 
clouds or photo-dissocitation regions (PDRs) at Ay < 2 (e.g. 
Keene et al. ll998l their Fig. 2). 

Ultra-high-resolution spectroscopy of CH + (Crawford et al. 
11994 1 absorption towards £Oph reveals broader profiles com- 
pared to CN and CH, with an upper limit kinetic temperature 
of ~2000 K, and with no velocity offset. This is consistent with 
CH + production through the endothermic reaction 

C + + H 2 + 0.4 eV CH + + H. (3) 

The analog of Eq. [5] shows that temperatures > 1000 K are 
required for CH + production via Eq. [5] CH + thus arises in 
the warm halo, or photo-dissociation region, surrounding the 
£Oph cloud. Crawford \\995\ and Crane et al. dl991t extend 
the £Oph results to five and twenty other lines of sight through 
translucent clouds (i.e. with Ay < 2). With such a kinetic tem- 
perature, it is clear that CO fractionation in the sites of CH + 
production is negligible. 

The carbon isotope ratio as measured by 12 CH + / 13 CH + is 
not affected by selective dissociation. Although photodissocia- 
tion of CH + contributes to PDR chemical networks (Sternberg 
& Dalgarno 1992), CH + is optically thin to photodissociat- 
ing UV radiation in translucent clouds. The strongest optical 
CH + line in <fOph is A l Tl - X l Y + at 4232 A (Morton 
which reaches maximum opacities of ~0.3 (e.g. this work). 
Since no ultra-violet CH + lines have been reported, they are 
bound to be faint, and thus thin. The photodissociation cross- 
section <t(v) calculated by Kirby et al. ( 1980) is mostly due to 
transitions from the ground state to the vibrational con- 
tinua of the 2'S + , 3 : S + , and 2'n states, which smooths out 
any difference between isotopes. Some resonant absorption de- 
rives from the calculation of cr(A) at wavelengths of 1418A, 
1453A, and 1490A. But Kirby et al. ( fT980l assign oscilla- 
tor strengths of less than 10~ 3 in any of these peaks, which 
represents Einstein B absorption rates one order of magnitude 
smaller than at 4232 A. 

The 12 CH + / 13 CH + ratio is thus a secure representation of 
12 C/ 13 C in translucent clouds. In addition, the absorption line 
from 12 CH + at 4232A is separated from the line of 13 CH + by 
0.265 A, which permits unblended measurements. These ratios, 
obtained from optical measurements, are restricted to regions 
within ~2 kpc of the Sun because of the need for bright back- 
ground stars, and of narrow velocity profiles in the intervening 
cloud. Eight ratios have been measured but only one line of 
sight, toward fOph, has an excellent signal to noise ratio. The 
average 12 C/ 13 C ratio from CH + data agrees well with those 
obtained for sources near the Sun, from radio astronomy data. 
This gives one confidence in the ratios obtained from radio 
data. 

However, the CH + measurements so far resulted in iso- 
topic ratios with a large scatter between different lines of sight, 
and there is only one line of sight with an excellent signal to 



noise ratio. The previous data from CH + were taken with the 
ESO CAT, and 4-m or smaller aperture telescopes. The best 
results from the point of view of the signal to noise ratio and 
the reliability of the ratio for £Oph were presented by 3 differ- 
ent groups (Stahl et al. [T989l Stahl & Wilson [T9921 Crane et 
al. rrWTl Hawkins et al. [T985l[TW3l . The final ratio obtained 
by each group was 67+5. The £Oph value is by far the best 
measurement of the 12 C/ 13 C ratio for the ISM near the Sun. 
However, the measurements for 7 other nearby weak sources 
were of lower quality. The ratios for 3 of these regions, to- 
ward HD26676 (64+6, Centuri on & Vladilo lTWTl . HD1 10432 
(71 + 11, Centurion et al. 119951 ). and //Norma (67 ±6, Centurion 
& Vladilo [l991l ) are in very good agreement with the £Oph 
ratio. But there are discordant ratios: R = 38 + 12 toward 
HD157038 (Hawki ns & M eyer 1989), and R = 49 + 15 toward 
£Per (Hawkins et al. [T993l while the ratio toward HD152235 is 
R = 126 ±29 and isR = 98± 19 towards HD152424 (Vladiloet 
al.[T993). Most of the measurements are limited by noise, but 
the profile toward HD 152424 is complex and thus a determi- 
nation of the ratio is not easy, because of CH + and 13 CH + line 
blending. 

We repeated the CH + observations toward £Oph, 
HD1 10432, HD152235, HD152424 and HD157038, and 
included the new lines of sight toward HD152236, HD154368, 
HD 16 1056, HD 169454, HD 170740. Our goal was first to 
determine whether the ratios are truly different from the value 
for ^Oph, which would indicate that the mixing of the ISM is 
not very fast and complete, and second, we wanted to increase 
the sample of lines of sight with good to excellent signal to 
noise ratios, and from this determine an average ratio near the 
Sun truly representative for the local ISM. 

Section |3 describes the observations, and Section [3] gives 
details on our procedure to measure the isotopic ratio R. 
Section[4]contains our results , which are then summarised and 
discussed in Section[5] Section[6]concludes. 

2. Observations and reduction 

There are two systems of CH + lines which have appreciable os- 
cillator strengths. Because of atmospheric absorption, the lines 
at 4232 A are easier to measure. The 13 CH + line is 0.26 A 
(Av = 19 km s~') shortward of the 12 CH + line. The other sys- 
tem, at 3957 A, is more difficult to measure, but has the prop- 
erty that the 13 CH + line is on the long wavelength side of the 
12 CH + line at a separation of 0.44 A(Av = -33 km s" 1 ). This 
property has been used by Stahl & Wilson ( 1992 ) to check that 
there is no accidental overlap of a weak velocity component of 
12 CH + at -18.8 km sec -1 from the deepest 12 CH + absorption 
line with the 13 CH + line. Nearly all previous data were taken at 
4232 A. 

Our observations were obtained with the UVES echelle 
spectrograph at the VLT unit telescope Kueyen at Cerro 
Paranal, Chile, in three nights between June 14/15 and June 
16/17, 2003. UVES allows measuring both line systems in the 
same detector setting, thus providing a means to correct for line 
blending if it is apparent in the spectra. 

The observations are difficult because of the need for both 
very high S/N and very high spectral resolution to acquire the 
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Fig. 1. Spectral region about the CH + absorption towards the 
template star HD205637, emphasizing the absence of telluric 
features at airmasses of 1.05-1.08. 

profile of faint 13 CH + , against a very bright stellar continuum 
and close to the much stronger 12 CH + component. We need 
S/N ratios of at least ~10 on 13 CH + for accurate profile fitting. 
The project requires the highest possible spectral resolution. 
We therefore used an image slicer to minimize flux losses. In 
addition, the image slicer distributes the light along the slit, 
which improves flat-fielding and allows for our bright targets 
longer integration times before saturation occurs. UVES slicer 
#2 was used, which reformats an entrance opening of 1'.'8 x 
2'.' 2 to a slit of 0744 width and a length of 7" 9, which is imaged 
on the spectrograph entrance slit of 0745 width. The spectral 
resolution in this configuration is A/AA = 75 000. The central 
wavelength was set to 4 370 A, which gives a spectral cover- 
age from 3 730 to 5 000 A. During our observations the seeing 
was typically too good to fill the entrance aperture of the im- 
ages slicer. Therefore, for most spectra, only one or two slices 
contained most of the signal, which unfortunately decreases the 
expected gain of the image slicer. 

Accurate flat-fielding is also important. Therefore a large 
number of flat-fields (200) was obtained during daytime dis- 
tributed along the observing run. 

In addition, a rapidly rotating unreddened early-type star 
(s Cap = HD205637, spectral type B3V) was observed in order 
to check for possible faint terrestrial features. We confirm there 
are no detectable telluric features under the CH + absorption, 
and show in Fig.^the template star spectrum in the region of 
interest, at airmasses < 1.08. Higher airmasses result in telluric 
absorption features, as in the case of £Oph discussed below. 

The brighter targets where exposed until about 50% of the 
maximum level allowed by the CCD detector was reached. For 
£ Oph, this limits the exposure time to about 5 sec, and typi- 
cally a few minutes for the fainter targets. Series of up to 50 
exposures per night were obtained to build up the required 
S/N-ratio. Some targets turned out to show too complex line 
profiles. These were dropped from our initial target list. The 
observations are summarized in Table ^ which also lists the 
mean airmass of the target objects in the different nights. At 
the higher airmass values, the range in airmass differed from 
the mean by about ±0.1 during the observations, and less at 
smaller airmass. 

We also list in Tabled the radial velocities relative to the 
solar barycentre of the CH + absorption. These are calculated 



from the shift in wavelength between the rest wavelengths of 
the transitions and of the average of the Gaussian centroids 
weighted by their equivalent widths. The uncertainty on the 
derived velocities is conservatively ~ 1km s _1 , and depends 
on the accuracy of the rest wavelengths. 

We used the Midas package originally developed for the 
ESO Feros spectrograph for the reduction of the spectra (Stahl 
et al.. U999> . In order to maximize the S/N-ratio of the extracted 
spectra, all flat-fields obtained during the run were averaged. 
After background subtraction and flat-fielding, the spectra were 
extracted with a very long slit, extracting all slices together. 
The wavelength calibration of the extracted spectra was done 
with a 2D-polynomial, fitting all echelle orders in one step. A 
mean ThAr-spectrum obtained during day time was used for 
the calibration of all spectra obtained in one night. Finally, all 
spectra were merged to a ID-spectrum and all observations of 
each night averaged to a nightly mean spectrum. Wavelengths 
are reported in air and refererred to the solar barycentre. 

Table 1. Summary of the observations, giving the total number 
of exposures, N, the integrated exposure time in seconds, and 
the observed heliocentric velocities of the absorbing clouds, 



object 




air mass 




N 


exp. 


v G 




14/06 


15/06 


16/06 




time 


km s _1 


<ro P h 


1.16 


1.81 


1.71 


150 


750 


-15.1 


HDl 10432 


1.31 


1.29 


1.29 


72 


6770 


+5.1 


HD 152235 


1.07 


1.31 


1.24 


65 


16250 


-4.3 


HDl 52236 




1.14 




10 


600 


-6.9 


HDl 52424 


1.05 






10 


2 500 


-7.6 


HDl 54368 


1.05 


1.06 


1.05 


50 


9 850 


-6.0 


HD157038 


1.11 






10 


2 500 


-14.7 


HD161056 


1.25 


1.06 


1.06 


59 


11 110 


-11.4 


HDl 69454 


1.30 


1.12 


1.12 


50 


22900 


-10.5 


HDl 70740 


1.05 


1.62 


1.60 


65 


7 800 


-11.0 


HD205637 


1.05 


1.08 


1.08 


150 


3 385 





The nightly spectra were combined in a weighted average 
to produce coadded spectra. The weights were taken as 1/cr 2 , 
where cr is the noise in each spectra, as calculated from the 
root-mean-square deviations from a linear fit to a region of the 
spectrum devoid of conspicuous features (we chose 4204. 2A to 
4205. 9A). Table|5]summarise the resulting weights. 

Table 2. Relative weights used to combined the nightly spectra 





June 14 


June 15 


June 16 


<ro P h 


0.45 


0.30 


0.26 


HDl 10432 


0.26 


0.37 


0.36 


HDl 52235 


0.41 


0.31 


0.28 


HDl 54368 


0.35 


0.39 


0.26 


HD161056 


0.31 


0.41 


0.29 


HDl 69454 


0.37 


0.34 


0.29 


HDl 70740 


0.31 


0.40 


0.29 


HD205367 


0.41 


0.32 


0.27 
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3. Model line profiles 

In order to use the fact that both 12 CH + and 13 CH + lines have 
the same opacity profile t(v), we must fit the 12 CH + absorption 
with a parametrised model and scale it to 13 CH + . 

We considered using Voigt profiles to account for the in- 
trinsic line profiles of the CH + vibronic lines. In order to fit the 
low-level broad wings (e.g. Stahl et al. fT9B9 1 towards £ Oph 
with an hypothetical Lorentzian core, Einstein A values of or- 
der 10 9 s _1 are required. Gredel. et al. (1993 1 quote an os- 
cillator strength of f a = 0.00545, or A = 2.258 10~ 4 , for 
12 CH + A4232 II - £(0,0) and/ = 0.00331, or A = 1.568 10~ 4 , 
for 12 CH + A3957 II-Z( 1 , 0). We therefore used simple Gaussian 
profiles to describe the lines. 

We preferred to fit the two overtones separately, and thus 
obtain independent measurements with which to assess the role 
of systematics. A simultaneous fit of both overtones could have 
helped constrain the opacity profiles, which in our separate fits 
do not always share the same Gaussian components. 

The fitting algorithm is as follows. 

1. Extract a 1.5 A spectrum F{ A) centred on the CH + line, Ai < 

A < A 2 . 

2. Define a baseline F C (A) with a 4 - 17 th order Legendre 
polynomial to the data, as in Sembach & Savage (1992i, 
F C (A) = aiPAJA - A l )l{A 2 - A,)}. 

3. Fit for the {a,}[ =1 coefficients in a least square sense, ignor- 
ing the neighbourhood of the ISM and telluric lines, and 
weighting the data to improve the quality of the baselines 
near the interstellar absorption lines. Store the rms disper- 
sion of the residuals as the spectrum's noise, cr f . o> is up- 
dated after fitting the parametrised model; it is replaced by 
the dispersion of the residuals in Step|8] 

4. Define a model spectrum with 

F*(A) = F c (A)cx V (-t(A)), (4) 

where the line absorption opacity t is a superposition of 
n g = 1 to 6 Gaussians on each isotope, with 13 CH + 
components sharing the parameters of the 12 CH + compo- 
nents, except their centroids are translated by a fixed ve- 
locity shift Av; so , and their amplitude are scaled by a factor 
/ = 13 CH + / 12 CH + , which is a free parameter in our fit: 

r(A) = J] [r° exp (-0.5(A - A°fl<r°]) + 

/ r° exp (-0.5(A - A°{\ + Av iso /c)) 2 /cr° 2 )] . (5) 

The value of Avj so depends on the overtone and is adjusted 
iteratively on £Oph and checked for consistency on all tar- 
gets. We use Av iso (3957) = 32.721 km s"\ Av iso (4232) = 
-18.967 km s -1 , with uncertainties of ~0.4km s _1 (roughly 
1/10 the resolution element). 

5. Convolve the model spectrum with the instrument re- 
sponse, B A : 

F m (A) = j dA'F*(A')B(A-A,). (6) 

We estimate the instrumental response from the arc lamp 
spectrum. The smallest line-widths are indicative of resolv- 
ing powers of « 83 000 ± 2000 rather than the nominal 



75 000 for our setup. The instrumental response is approx- 
imated to a Gaussian, with a FWHM given by /l o /83000, 
where A a is the center wavelength of each CH + overtone. 

6. Calculate the goodness of fit function 

X 2 = J](F(Aj) - F m (Aj)) 2 /cr 2 F , (7) 

using the noise from Step [3] We use a single noise value 
and neglect possible variations of the noise level, as would 
be the case for Poisson statistics. The noise variations are 
important for very steep spectra, which is not our case, or 
under deep absorption lines. But the accuracy of the model 
line profiles under 12 CH + is generally not important to infer 
isotope ratios (only the wings matter). 

7. Optimize the model parameters by minimizing^ 2 . Perform 
an initial estimate of the global minimum with a genetic 
algorithm (pikaia, Charbonneau 1995 1, and optimize with 
the downhill simplex method (amoeba, Press et al. 1986). 

8. Update the noise of the stellar spectrum by replacing cr F 
with the rms dispersion of the residuals (the difference be- 
tween the observed and model spectra). 

9. Estimate the uncertainty in model parameters by two meth- 
ods: 

(a) search parameter space one parameter at a time to find 
the region enclosing 68.3% confidence for the^ 2 distri- 
bution with one degree of freedom, thus obtaining dif- 
ferent upwards and downwards error bars and account- 
ing for statistical bias, 

(b) estimate lcr uncertainties from the curvature matrix of 
X 2 (i.e. approximate to normal errors). 

The UVES spectra are oversampled, so that the spectral dat- 
apoints are correlated. Independent datapoints are separated by 
roughly one resolution element, of AA = /1/75000 ~ 0.05 A. 
Since the spectra are sampled with an interval of 0.01 A, 
Ncor ~ 5 consecutive datapoints are correlated. In some cases 
the number of Gaussian components involves as many free- 
parameters as independent data points (as for HD 154368 and 
HD169454). Some Gaussian components may thus be redun- 
dant. But we caution that our purpose is to obtain isotopic ra- 
tios, not to report on the velocity structure of the ISM in CH + . 
To constrain the opacity profile consistently we would have had 
to fit the two overtones simultaneously. 

The correlation between the spectral datapoints is ignored 
when fitting for the parametrised model. The off-diagonal terms 
of the covariance matrix add terms to the x 2 goodness-of-fit 
estimator (Eq. 0. The number of additional terms is roughly 
N CO i per diagonal element, which increase^- 2 by a factor N cor . 
But the noise estimate from the dispersion of the residuals also 
increases due to the correlation of the datapoints, the actual 
noise should be o> x ~\/N cor . Both corrections cancel out and 
leave x 2 as in Eq.[7] 

We tested for systematic noise by comparing the noise in 
the nightly coadded spectra, cr co ^d, with that expected from 
the noise in each of the N exp individual exposures, en . For ther- 
mal noise it should hold that <x coadd - cr\l sjN exp . We define 

Q = o"i /( -yWexpCcoadd), an ^ ra i se tne 01 "der of the baselines un- 
til Q ~ 1. The tests were run for the case of HD 170740, which 
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Fig. 2. The ratio of measured to expected noise, Q, as a func- 
tion of the order of the Legendre polynomial used to define 
the baselines, /. The data are for 25 exposures of HD 170740 
on June 15 th , for two 1.5 A spectral regions centred on AA232 
(solid line), and on /13957 (dotted line). 



presents the most ragged continuum. The reference baseline 
is defined on the coadded spectrum, and the noise is given 
by the dispersion of the residuals, excluding the spectral re- 
gions affected by the ISM lines (these regions are defined on 
Fig. I12> . We scale the reference baseline to the line-free con- 
tinua of single exposures using the ratio of their median values. 
Fig.|2]shows we need to reach orders of / ~15 for A3951, and of 
I ~8 for AA232. Below these orders the continua are affected by 
systematic noise. Unfortunately the level of systematic noise 
changes for each observation and setting, so that we can only 
use the tests on HD 170740 as guidelines for the other lines of 
sight. In practice we chose the smallest order that is compatible 
with the data, without exceeding the limits in HD 170740. 

The isotopic ratio derived from the fits is the ratio of the 
column densities of each isotope, R = N( l2 CU + ) /N( l3 CU + ). 
With the best fit / we can calculate R = l/f for a line of sight 
with constant R. We neglect the small difference in oscillator 
strength between 12 CH + and 13 CH + . The two ions have slightly 
different vibrational structures which will lead to small changes 
in the /„ values for individual vibrationally resolved transitions 
(J. Tennyson, private communication). We are thus neglecting 
f -value differences of the order of the relative difference in 
reduced mass, or A/ // ~ 1 /170. 

In taking statistics on the values ofR - 1 // in the local ISM 
care must be taken to assess the significance of measurements 
along individual sightlines. The uncertainties on individual fits 
do not include the systematic error involved in baseline defi- 
nition. If baselines were known a priori, the expectation value 
and uncertainty on / for each individual sightline would be ob- 
tained from the statistics of multiple measurements of the same 
quantity. Thus, in the absence of this systematic error, (/) and 
<r(f) correspond to the weighted average and quadratic sum of 
the weighted uncertainties, where the weights are taken as the 
inverse variance of each measurement, 1/cr 2 : 

N 

</> = Yjf m > (8) 



N N 

°!(f) = Tj {(Ti wi)1 = vYp-iTi?* ( 9 ) 
;=i i=i 

for N measurements on a given sightline, and where the indi- 
vidual errors are those derived from the curvature matrix. 

To estimate the systematic error involved in baseline defi- 
nition, we compare a posteriori the scatter of individual mea- 
surements with that expected from the noise level. We use the 
rms dispersion 

<r 2 (f) = [(/,-</>)] 2 . (10) 

The uncertainty o~\{f) neglects the systematic errors. <x 2 (/) 
uncertainties will give larger error bars, but are probably too 
conservative considering we are in fact repeatedly measuring 
the same quantity. Thus the actual uncertainty on our measure- 
ment of / for each sightline is intermediate between cr\ (/) and 
o~i(f)- For this reason we list both values. We have set <x 2 = <J\ 
if <x 2 < <T\, as is the case when only one ratio is available to 
sample the effects of systematic uncertainties in the continuum 
definition (or when one ratio has much higher weight than the 
others). 

4. Results 

An application of the fitting procedure described in Sec. [3] to 
our sample of 10 lines of sight gives the best fit parameters and 
isotopic ratios listed in Table[3] 

The isotopic ratio / = 13 C/ 12 C is better suited as free pa- 
rameter than R = 1 // for the purpose of fitting individual spec- 
tra, and subsequently averaging the best fit values. For noisy 
data the error propagation when using R would involve second 
order expansions in R(f) when calculating averages and disper- 
sions. In what follows we prefer to list our results in terms of /, 
and revert to the more common usage of R in Section [5] when 
comparing with previous data. 

Table lists the isotope ratios and best fit parameters. It 
is apparent the AA232 line generally gives the best estimate of 
R. We note from Table [3] that the normal errors are good ap- 
proximations to the 68.3% confidence limits, which justifies 
their subsequent use in the weighted averages. We also note the 
weighted average of several measurements for the same line of 
sight are consistent, within the errors, with the fit to the coad- 
ded spectra (which are in general excluded from the averages), 
giving confidence in our procedure. 

4.1. Notes on individual objects 
4.1.1. fOph 

From Hipparcos results, the distance to £Oph gives an upper 
limit to the absorbing cloud of 140±14 pc. The spectra of £Oph 
around the CH + lines are shown on Fig. [5] where we have also 
plotted the Gaussian components on the 12 CH + line (the indi- 
vidual components on the 13 CH + line are left out for clarity). 

Residual features under the absorption lines vary on scales 
of ~0.02A . These reflect the imperfection in the model line 
profiles, more Gaussian components reduce these residuals. 
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Table 3. We list / x 10000 to improve the clarity of the table. cr(f) refers to the Gaussian errors (those derived from the 
X 2 curvature matrix). The average values are quoted as (/) ± cr 1 (cr 2 ) (see text for details), and exclude the coadded spectra. 
Equivalent widths for 12 CH + , W A , are given in mA. The stars select the measurements included in the combined values of R; the 
selection are discussed in the text. 
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The £Oph spectra were acquired at airmasses of 1.8 and 1.7 
on the 15 th and 16 th of June. Both spectra show absorption at 
3958. 15A (at 3958.5 A in the observatory's rest frame). This is 
absent in the spectrum from 14 th of June and in all other tar- 
gets. Since all other spectra were acquired with low airmass 
values compared with that of £Oph on the 15 th and 16 th , we 
conclude the 3958. 15A feature is telluric. We did not attempt 
to include the 3958. 15A feature in the fits. The telluric absorp- 
tion propagates into the residuals, and contributes to the noise 
level. Note that the telluric feature has no effect on where we 
set the baseline level, because we assigned zero weights to the 
spectral points in its neighbourhood. 

Stahl et al. c l989J reported R = 77 ± 3 1 from ,14232, us- 
ing a single Gaussian fit. But Stahl & Wilson ( 1992 1 used two 
Gaussian components to account for the broad wings detected 
by Crane et al. ([1991 1, and thus obtained a tighter fit to the 
^OphCFL absorption, wifhfl = 71 ±3 for ,14232 and/? = 68±7 
for ,13957. 

We also find that using two instead of one Gaussian com- 
ponents significantly improves the fit, bringing reduced^- 2 val- 
ues for the coadded £Oph ,14232 spectrum from 2.07 for one 
Gaussian (using the noise reference of the two Gaussian fit), 
to 1.05 for two Gaussians. But the resulting average for R 
is unchanged, we obtain (R) = 78.22 ± 2.63(3.75) with one 
Gaussian, {R) = 80.20+2.75(4.19) with two Gaussians. 

What is the effect of including more than two Gaussians? 
It is difficult to separate an error in the definition of the con- 
tinuum baseline from a real increase in confidence level of the 
fit. We experimented on £Oph with 3 and 4 Gaussians, keep- 
ing the noise fixed to the rms dispersion of the residuals in 
the 2 Gaussian fit. The reduced^ 2 were respectively 0.92 and 
0.88 for the 3 and 4 Gaussian fits to the spectral region cen- 
tred on /14232 in the coadded spectrum of £Oph. But some 
of the 4-Gaussians fits to the ( Oph spectra contained very 
narrow features, corresponding to unphysical opacity compo- 
nents. Since we cannot separate the uncertainty in the baseline 
definition, we adopt the smallest number of Gaussians consis- 
tent with the dataset (with reduced x 1 close to 1). We chose 
3 Gaussians rather than 2 because Crawford et al. 1994 report 
a 3-Gaussians decomposition of the ultra-high resolution CH + 
profile toward £Oph. The final value from the UVES data is 
thus (R) = 78.47 ± 2.65(3.53), with a 3-Gaussian fit. 

All R values for £Oph are consistent within 2<x. The value 
we obtain here is higher than previously reported by Stahl et al. 
( 1992 1 by 1.6 cr, adding the statistical errors in quadrature. For 
comparison Hawkins et al. ( 1993 1 give R = 63 ± 8 for /14232, 
67 + 19 for /13957, Crane et al. fTWil give R = 67.6 + 4.5 for 
,14232, and Vanden Bout & Snell JT9801 report R = 77+ \ 7 2 , also 
for ,14232. 

The value reported in this article has an accuracy close to 
that of Stahl et al. ( 1992|l, of 71 + 3. Combining both measure- 
ment with identical weights gives our best value for £Oph as 
R = 74.7 ± 2.3. 

In what follows we will nonetheless use the UVES value 
of R for £Oph to compare with the other lines of sight, (R) = 

1 note their uncertainty is estimated as for o"i in this work, in the 
case where all samples are assigned the same weight 



78.47 ± 2.65(3.53) thus relying exclusively on measurements 
obtained with the same instrument. 

4.1.2. HD1 10432 

The profile from Crawford ( 1995 1 is well fit by 1 Gaussian only, 
but we obtain best results with 2 Gaussians. 

The Hipparcos distance to HD1 10432 is 300±51 pc, which 
confirms that the CH + absorption most likely arises in the 
Coalsack, at a maximum distance of 180 pc (Franco 1989 1. 

The fit on the ,13957 line is excluded from the combined 
measurement of R because the adjacent telluric absorption fea- 
ture could affect the wings of the CH + absorption. Also the 
baseline around ,13957 seems to be strongly affected by sys- 
tematic noise. The ,13957 values differ from ,14232, which oth- 
erwise give consistent numbers for each night. 

4.1.3. HD1 52235 

HD 152235 is a member of the Sco OBI association, but as 
discussed by Crawford dl995> the absorbing CH + lies in the 
Lupus molecular cloud, at 170 pc (Murphy et al. 1986 1. 

The profile from Crawford ( 1995 1 is fit by 3 Gaussians. Our 
data, presented in Fig.|5] also required 3 Gaussians for ,13957. 

The CH + absorption has a red tail than difficults separating 
both isotopes at 4232A, while the telluric absorption feature 
borders the 13 CH + line. It is difficult to pin down the under- 
lying continuum at 4232A. We tried two approaches to define 
the continuum from edge to edge of the CH + absorption (from 
4232A to 4232. 7A): a free fit (i.e. roughly a straight-line, with 
zero weights over 4232A - 4232. 7A), and the inclusion of a 
small non-zero weight mid-way between the two isotopes (see 
Fig.0. 

The spectrum from 2003 June 14th is clear of telluric 
absorption (it was observed with the smallest air mass, see 
Table [0, so both overtones should give consistent values for 
R. This is why we chose the baseline shown on Fig.|5] The two 
overtones give consistently high value for June 14th, and also 
on average. So we are rather confident of the ratio we report. 
The fits on ,13957 for June 15 th and 16 th are excluded from the 
combined measurement of R because they are affected by tel- 
luric absorption under the red tail of I3 CH + . 

Vladilo et al. ( 1993 1 reported R = 126 + 29. Our measure- 
ment, which has better accuracy, gives lower value, although 
still higher than the ISM average (see Sec.|5}. 

This line of sight gives a significantly higher R value than 
in HD1 10432 at 4.2 cr, using the conservative uncertainties 0"2, 
from the scatter of each individual measurement of R. 

AAA. HD1 52236 

Fig.[6]presents the CH + spectra for HD152236, along with our 
fits. Crawford 0995 1 used 4 Gaussians, but we needed 6. 

HD152236, as HD152235, is a member of the Sco OBI 
association, so that the absorption probably occurs in the Lupus 
cloud. 
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Fig. 3. £Oph spectra. Wavelengths are in Angstroms, and flux densities in arbitrary units. The individual Gaussian components 
comprising the fit, prior to folding with the instrumental response, are shown only on 12 CH + as light gray lines, but are omitted 
from 13 CH + for clarity (see text for details). The vertical dashed lines mark the line centroids, at the average of the Gaussian 
centroids weighted by their equivalent widths. The units of the y-axis are arbitrary, and are scaled so that the median of the object 
spectrum is unity. Also shown is a magnified version of the object spectrum, by the factor indicated on the figure, and offset for 
clarity. The residual spectrum is shown on top of the object spectrum, and is also magnified by the factor indicated on figure. 
The noise used to compute the significance of the fits is labelled on the top right. The noise is the rms dispersion of the residual 
spectrum, and is in the same units as the j-axis. The height of the shaded rectangles on top of the x-axis indicates the relative 
weights used in the baseline definition. 
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Fig. 4. Same as Fig. [3] The fits on the A3957 lines are excluded from the combined measurement of R. 

There is no manifest absorption at A3957, but we give the Although the A4232 line has better signal to noise, both 

result of the formal fit, which allows assigning a lower limit R isotopes are blended. We tried to apply the technique described 

value. We take /upiimit = (/) + 3cr 1 (/), and/? > l// up iimit, giving below in the case of HP 154368 (Sec. l4.1.6> . but did not obtain 
R > 61.4 at the 3 cr level. 
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Fig. 5. Same as Fig. [3] The fits on A3951 for June 15 th and 16 th are excluded from the combined measurement of R. 



significant excess absorption in the blue wing of the main line. 4.1 .5. HD1 52424 
The lower limit refers only to the /13957 measurement. 



HD152424, as HD152235, is a member of the Sco OBI as- 
sociation, so that the absorption probably occurs in the Lupus 
cloud. 
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Vladilo et al.Q993 report/? = 98 ± 19, but this target gives 
us our lowest R value, of 55.65 ± 3.55(15.61). 

The two isotopes are blended at /14232, so we apply the 
same technique as for HD154368 (Sec. I4.1.6> . The residu- 
als under the 13 CH+ absorption are consistent with the noise, 
which gives us confidence in the inferred R value. HD 152424 
provides our best measurement of R for blended lines, as ob- 
tained by keeping the opacity profile fixed to the best fit on 
-13957. But defining consistent baselines between the two over- 
tones is very difficult without including the baseline as free pa- 
rameter. So we decided to ignore the significant residuals under 
the main isotope at -14232. We include the noisy residuals in the 
noise estimate, so that the significance of the R values derived 
from the blended line is lowered. 

The importance of quantifying systematic uncertainties is 
manifest in HD152424: a visual inspection of Fig. shows 
the fit under A3951 is very good, and that the baseline seems 
smooth. However appeareances can be deceiving: the extrap- 
olation of the -13957 fit to /14232 gives a very different value 
of R. Had we considered only the cr\ uncertainties, the value 
for HD152424 would have been one of the most accurate. But 
because CH + towards HD 152424 is so broad the conservative 
uncertainty (j% is much higher than cr\ and reflects the system- 
atic uncertainty involved in the baseline definition. 

4.1.6. HD1 54368 

This target was selected from the survey of Gredel et al. ( 1993 1. 

The CH + absoprtion for HD 154368 (with an Hipparcos dis- 
tance of 364.8 + 128.3 pc) shown on Fig.[S]has four conspicu- 
ously distinct velocity components. The two isotopes are man- 
ifestly blended for the /14232 transition. 

Our data cover the two overtones in the same UVES 
dichroic setting simultaneously, and the spectral resolution is 
expected to be constant with wavelength in such an instrumen- 
tal setup. After checking we had no detectable trend of varying 
resolution by measuring the arc line widths, we attempted fit- 
ting /14232 using the results of the /13957 fits. 

A fit of both overtones simultaneously would be counter- 
productive because the -14232 region has better sensitivity, and 
would dominate the fit, placing Gaussian components of 12 CH + 



under the rarer isotope. Instead we kept the 12 CH + opacity pro- 
file as a function of velocity r(v) as inferred from the /13957 
fit. We then scaled t(v) to t(-1) for the -14232 region, which 
implies scaling component widths and separations in wave- 
lengths. After setting baselines as in the unblended case, the 
optimization involved two free parameters, the scaling factors 
on the 12 CH + and 13 CH + opacity profiles. 

The fits shown on Fig. [3] are rather poor in the case of 
-14232. Uncertainties in the baseline determination resulted in 
significant residuals under the main isotope line. But these do 
not affect the 13 CH + absorption. We recompute the noise level 
from the residuals, including the baseline uncertainties. The un- 
certainties on R tabulated in Table [5] are a result of this exag- 
gerated noise level. 

4.1.7. HD1 57038 

This line of sight also has blended absorption at /14232, and we 
apply the same technique as for HD 154368. But the results on 
-14232 are useless because they are too sensitive on the baseline 
definition, so we avoided the use of -14232 altogether. 

Hawkins & Meyer ( 1989 and references therein) give R = 
38 + 12 from /14232 and use a distance of 1 .7 kpc to HD157038. 
We confirm this line of sight is enriched in 13 C, with a lower R 
value than the average ISM value (see Section l5~H . 

4.1.8. HD161056 

13 CH + absorption towards HD161056 (selected from Gredel 
et al. I1993I is very weak at -13957, and yet conspicuous at 
-14232 (see Fig. llOb . We first attempted to fit the two overtones 
separately, and found widely discrepant values: /(3957) = 
68.9±11.8(12.0),and/(4232) = 159.3+5.6(19.1).The two iso- 
topes may be blended at /14232. We thus applied the same tech- 
nique as for HD 154368, extrapolating the -13957 opacity profile 
to -14232. But the inferred isotope ratios were still very different 
(see Tabled There must be additional absorption bridging the 
two overtones at -14232, not due to CH + : the red tail of -14232 
is also seen at -13957, while the absorption at the blue edge of 
-14232 is not seen at -13957. We cannot identify the nature of 
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Fig. 7. Same as Fig.[S] 13 CH + is visible as excess absorption on the blue edge of AA232 I2 CH + 
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Fig. 8. Same as Fig.|3] This line of sight has blended I2 CH + and 13 CH + absorption at AA232. We extrapolate the 12 CH + profile at 
/13957 (see text for details). 13 CH + is visible as excess absorption on the blue edge of I2 CH + . 
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Fig. 9. Same as Fig. [8] The fit for CH + absorption on the blue edge of /14232 CH + gives a 3 cr detection, but is barely visible 
on the Figure. 



the additional absorption, so we choose to exclude AA232 from 4.1 .9. HD1 69454 
the combined measurement of R in HD161056. 

This line of sight, selected from Gredel et al. d!993> . exibits 
at least 4 distinct velocity components. We use 6 Gaussians to 
describe the profile, but did not attempt to compare reduced^ 2 
when varying the number of Gaussians because of the uncer- 
The hipparcos distance to HD161056 is 425.7 ± 145.5 pc. tainties in the baseline definition. 
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The total number of free-parameters used to describe the 
main line is 6x3, with 3 parameters per Gaussian. The zero- 
absorption line width for HD169454 is -0.8 A, or 80 spectral 
datapoints. But the number of independent resolution elements 
under 12 CH + is about 16. Thus 18 free parameters may seem 
slightly excessive. But over-constraining the opacity profile is 
of no consequence to the derived isotopic ratio. 

We use the same technique as for HD 154368 to model the 
blended 12 CH + and 13 CH + absorption at A4232, but do not ob- 
tain satisfactory results. The fits on /14232 are shown on Fig. II II 
for completeness, and are not included in the combined mea- 
surement of R towards HD 169454. An unphysical value of / 
is obtained, but this value is still consistent with zero 13 CH + 
within the uncertainties. Also excluded from the combined 
measurement of R is the fit to /13957 from June 14 th because 
it shows non-ISM absorption at 3957.3 A. The ,13957.3 feature 
is not seen in other spectra acquired at similar airmass (e.g. 
that of HD 152235 on June 15 th ), and is either an instrumental 
artifact or the effect of passing clouds. 

4.1.10. HD1 70740 

This line of sight, Selected from Gredel et al. ( 1993 ), is impor- 
tant since it has the highest R value among those considered in 
this work, together with HD161056. The Hipparcos distance to 
HD 170740 is 21 1.9 + 42.4 pc. 

Unfortunately instrumental fringing, especially at ,13957, 
limits the accuracy of the fits. The ,13957 I3 CH + line is only 
marginally detected in individual nights, but we give the formal 
fits. The value we report relies on the nigthly spectra for ,14232, 
and on the coadded spectrum for ,13957. 

5. Discussion 

5.1. Observed scatter and comparison with other 
sightlines 

Figure[T3lsumm arises our results. The isotopic ratio R is lower 
for HD1 10432 than for HD152235 by 4.3 cr, and lower than 
for HD 170740 by 3.2 cr. In this comparison we have used the 
conservative uncertainties cr 2 , derived from the scatter of R for 
different nights and different overtones. 

The weighted average of R for the 9 lines of sight consid- 
ered in this work with 13 CH + detections is Rism = 12 C/ I3 C = 
78.27 ± 1.83. The weighted 1 cr dispersion of the data is 12.7. 
The hypothesis that the observed scatter is derived from a sin- 
gle value of R can be tested by calculating the corresponding 
value of ^ 2 /v, where v = 8 is the number of the degrees of 
freedom (the number of line of sights less one free parameter). 
We obtain x 2 / v - 6-0, which discards the hypothesis at a con- 
fidence level of essentially 1 . 

The observed scatter is not due to the galactocentric abun- 
dance gradient studied by Hawkins & Meyer ( 1989 1, and pre- 
dicted by Galactic chemical evolution models. The targets in 
our sample all lie within 30 deg of / = 0, and at distances of 
at most ~200 pc (except for the clouds toward HD157038 and 
HD169454, whose distances are unknown). Fig.[0]also gives 
longitude and distance information for each target. 
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Fig. 13. A plot of individual R values for the lines of sight stud- 
ied in this work. The thick lines are +cn uncertainties, from 
the weighted average uncertainty of individual measurements. 
The thin lines are ±0"2 uncertainties, from the scatter of mea- 
surements for a given line of sight, and are set equal to a"i if 
cr 2 < <T\. The hatched box is the weighted average of the mea- 
surements in the figure, and should be representative of the lo- 
cal ISM. This value is 78.27 + 1.83. The height of the hatched 
rectangle is ± 1 cr. Galactic longitude and distance in parsec are 
indicated to the bottom right of each target. 

Our value for R\su is significantly higher than those for the 
Pleiades, £Per and P Cyg. Vanden Bout & Snell fT9"8"0T l mea- 
sured 49+f and 59^ for 20 Tau (Maia in the Pleiades) and 
£Per. Hawkins & Jura ( 1987 1 report R values of 46+15 for £Per, 
40+9 for 20 Tau, 41+9 for 23 Tau (also in the Pleiades), and 
45±10 for P Cyg. Hawkins et al. ( 1993 ) confirm a low value of 
49±15 for £Per. Even taking 3 cr errors, the measurements by 
Hawkins et al. are below our results. One can ask whether the 
13 C enrichment in the Pleiades is due to contamination from 
past stellar winds in the cluster itself. 

Aside from the R values obtained by other groups and 
quoted so far in this work, there are also those from Centurion 
& Vladilo JT99H for //Nor (67+6) and HD26676 (64±6), 
which differ at 2 cr with our average ISM value (including its 
scatter). 

5.2. Are there D/H variations in the local ISM? 

From our data set, the values of R vary by 7 cr in the local 
translucent clouds, where cr is the uncertainty in the weighted 
average value (Rism)- This is strong indication that the ISM 
at the same galactocentric distance is not completely mixed. 
If this is the case in the chemical composition of the ISM 
in general, then one expects that the D/H ratios in the local 
ISM will also differ from place to place. Moos et al. (2002i 
(see also Vidal-Madjar 2002 1 summarise the FUSE results to- 
wards nearby white dwarfs, located within 100 pc of the Sun. 
The D i/H i ratios measured with interstellar D i and H i Lyman 
absorption varies by less than ~10% towards 7 white dwarfs 
within 100 pc of the Sun, at the limits of the Local Bubble. The 
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fractional standard deviation of D i/H i (the ratio of the scat- 
ter in the data to its weighted average) is 12%, but cannot be 
distinguished from a single D i/H i value within the uncertain- 
ties. By contrast, the D i/H i ratios obtained with IMAPS show 
a factor of 2 scatter towards y 1 Vel, £ Pup and 5 Oh, beyond the 
Local Bubble, at distances of 300-500 pc (Jenkins et al. 1999 
Sonnebornetal. l2000b . 

Is the scatter in the observed interstellar D i/H i due to vari- 
ations in D/H? The primary objective of the IMAPS Di/Hi 
observations, as stated by Sonneborn et al. (2000), is "to de- 
termine D/H with sufficient accuracy to test for spatial inho- 
mogeneity". But the D/H ratio is notoriously difficult to mea- 
sure, and often involves the use of several different instruments. 
The Hi lines are > 100 km s _1 broad and saturated so that 
N(H i) is measured from the Ly a damping wings, while N(D i) 
is best measured from Ly 6 or Ly e in the case of the IMAPS tar- 
gets. The D/H ratio involves carefully considering the effect of 
the underlying stellar continua and the contribution from other 
interstellar Hi components in the rising part of the curve of 
growth (with a high ratio of equivalent width to opacity com- 
pared to the main absorption line). Chemical fractionation and 
selective dissociation of HD affects D i/H i for lines of sight 
where H2 is observed, and differential acceleration by radia- 
tion pressure could displace Di relative to Hi in the diffuse 
ISM. Bruston et al. J198U assign the D i/H 1 scatter known at 
the time to such effects, and estimate a value of unperturbed 
D/H almost a factor of two higher than the observed average 
D/H of 1.5 10~ 5 . 

In addition Draine ( 2004 ) argues for variable depletion of D 
onto dust grains, such that the gas phase D/H would be a func- 
tion of dust processing, which is in turn a function of Galactic 
environment. The variations in D i/H 1 do not seem to be a faith- 
ful measure of D/H. 

5.3. Inefficient mixing in the local ISM 

What is the level of elemental homogeneity in the ISM? What 
is the power-spectrum of the dispersion in elemental abundance 
relative to the total metalicity? Does heterogeneity increase 
with size? 

At least one extreme of spatial scale indicates perfect mix- 
ing: the meteoritic evidence favours mixing at the molecular 
level of interstellar dust from different origin, during the for- 
mation of the solar system. Zinner et al. (119911 their Table 1) 
measured '2C/ I3 C values from 3 to 1000 in different SiC grains 
found in the Murchison meteorite. At some point in the presolar 
cloud the SiC grains must have been mixed from a previously 
heterogenous distribution. 

In a recent review of interstellar turbulence and mixing, 
Scalo & Elmegreen (2004i summarise current knowledge of 
the level of elemental heterogeneity in the ISM. The stellar 
return to the ISM should be spotty and poorly mixed. Scalo 
& Elmegreen explain that in contrast with diffusive processes, 
turbulence transport does not homogenize the gas at the atomic 
level. Yet the heterogeneity is notoriously difficult to detect. 
The available gas-phase diagnostics depend on local excitation, 



while the stellar data give upper limits only on the dispersion 
of elemental composition. 

But the 12 CH + / 13 CH + ratio, in contrast to Di/Hi, is a 
particularly sensitive probe of the mixing efficiency of the 
ISM. In terms of the fractional standard deviation used 
by Moos et al. (2002i, the scatter we have measured is 
o"ism/(^ism) =16.2% on a scale of 100 pc. Considering the 
precise location of the absorbing material is unknown, we take 
the Moos et al. ( 2002 ) measurement of a 10% scatter of D i/H 1 
from the mean value as a indication of inhomogeneity of the 
local ISM, over 100 pc scales. This value is similar to the vari- 
ations reported in this work, even though D is burned in stellar 
interiors, rather than produced as is the case of 13 C. Although 
their absolute yields are different, the mixing process (i.e. pas- 
sive scalar turbulence, Scalo & Elmegreen 2004) should be 
similar for both species. 

5.4. Is the solar value for R comparable to ISM 
measurements? 

The terrestrial value of R is 90 (Rosman & Taylor 1998), and 
is usually extrapolated to the Sun. However, the preliminary 
results of Ayres et al. (120051 ) indicate a photospheric carbon 
isotopic ratio of 70. The average value of our measurements of 
R in the local ISM, 78.27 ± 1.83, is slightly lower than the solar 
value of 90 (or even higher than the photospheric value of 70). 
But in fact (Rism) is indistinguishable from the solar carbon 
isotope ratio given the observed scatter of 12.7. 

Models of galactic chemical evolution predict an enhance- 
ment of 13 C relative to 12 C with time. Thus our result that 
13 C/ 12 C has essentially remained constant over the past 4.5 Gyr 
is surprising. The inconsistency is worsened when comparing 
with the preliminary results of Ayres et al. (2005). We may be 
affected by a sampling bias. Around 30 lines of sight or more 
are required to faithfully estimate (Rism)- 

But is the solar value truly representative of the ISM at the 
time of collapse of the solar nebula? The process of CO frac- 
tionation followed by condensation on cold dust grains may 
have increased 13 C in the pre-solar dust. With an enhanced 
dust-to-gas ratio through sedimentation in the accretion disk 
of the Sun, the result may have been a modification in the final 
R value. 

We imagine a 13 C-rich disk with a dust mass of 0.01 M , 
dust to gas ratio of unity, accreting on a zero-age Sun with met- 
alicity Z=0.01. The 13 C-poor atmosphere of the disk may have 
been blown away by the early solar wind, so that once the ac- 
cretion of the disk is concluded, the Sun reaches Z = 0.02 and 
higher 13 C/ 12 C than the pre-solar cloud. 

Aside from the possibility of preferential accretion of 13 C 
over 12 C, there is evidence that the solar nebula differed in 
composition from the bulk ISM. For instance the metalicity 
of B stars in Orion is lower than that of the Sun, even after 
4.5 Gyr of Galactic chemical evolution, indicating that the Sun 
must have been exceptionally metal-rich. This discrepancy has 
been interpreted in the framework of Galactic diffusion ( Wielen 
119771 ). by which the Sun diffused from its birthplace outwards 
in Galactocentric radius by ArcAL = 1-9 kpc. Since 13 C/ 12 C 



S. Casassus et al.: Interstellar 12 C/ 13 C 



15 



is also predicted to increase with tgal. Galactic diffusion may 
bring in agreement the ISM and solar values of R at the time of 
solar birth (Wielen & Wilson [TW7b . 

But perhaps the most compelling evidence that isotopic ra- 
tios in the Sun cannot be taken as direct constraint on the bulk 
ISM are isotopic anomalies in meteorites, such as the presence 
of extinct radionuclides . For instance Zinner et al . ( 1 99 1 1 found 
that 26 A1/ 27 A1»1 in the pre-solar nebula. They measured the 
concentration of the decay product of 26 Al, 26 Mg, locked in 
SiC grains from the Murchison meteorite. But in other mete- 
orites 26 Al/ 27 Ala; 1/4000 in Al-rich minerals from refractory in- 
clusions not representative of the pre-solar composition. Since 
the half-life of 26 Al is 10 6 yr, the pre-solar enrichment in 26 Al 
suggests pollution by a nearby source of 26 Al, such as a super- 
nova explosion (see Zinner et al. 119911 and Clayton 119941 for 
thorough discussions). 

6. Conclusions 

The VLT-UVES spectra of CH + absorption in the lines of sight 
toward 10 bright stars have allowed us to measure the carbon 
isotopic ratio R = 12 C/ 13 C with unprecedented accuracy. We 
confirm previous measurements of R toward £Oph, and obtain 
significant scatter in the local ISM. 

Our value for ^Oph is R = 78.47 + 2.65(3.53), where the 
second uncertainty in parenthesis is that obtained a-posteriori, 
from independent measurements, and independent baseline and 
absorption fits. Combining this value with that of Stahl et al. 
JT9921 gives our best value for £Oph, R = 74.7 + 2.3. 

Averaging our measurements for the 10 lines of sight gives 
a value representative of the local ISM: (/?ism) =78.27 + 1.83, 
with a weighted rms dispersion of 12.7. The dispersion in iso- 
topic ratio is 7 times the uncertainty in (Rism) - we detect het- 
erogeneity at 7 cr. 

The observed scatter in 13 C/ 12 C is the first significant de- 
tection of heterogeneity in the isotopic composition of the local 
ISM, and was obtained with a single instrument and a homoge- 
neous analysis, with a diagnostic independent of the local ex- 
citation and of chemical fractionation. The observed variations 
in R can be extrapolated to an overall elemental heterogeneity, 
and reveal that mixing in the ISM is not perfect. 

The solar carbon isotope ratio of 90 is undistinguishable 
from the present-day ISM ratio, considering its intrinsic scat- 
ter, even after 4.5 Gyr of galactic evolution. This statement is 
not affected by instrumental uncertainties, although a sampling 
bias could distort our ISM value. We question whether the solar 
carbon isotopic ratio is equal to that of the presolar nebula. 

A larger sample of stars, and new data towards the Pleiades, 
the Taurus cloud and P Cyg with 8 m aperture telescopes, will 
allow a better estimate of the average ISM value, and an im- 
proved understanding of the relationship between elemental 
heterogeneity and spatial scale in the ISM. 
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Fig. 10. Same as Fig. [3] None of the A4232 spectra were used to combine individual measurements of R. There is excess ab- 
sorption bridging the two isotopes at ,14232. The origin of the additional absorption, whether a telluric or a stellar feature, is 
unclear. 
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Fig. 11. Same as Fig. [8] The information from A4232 is not used because / is spuriously negative, and neither is that from /13957 
on June 14 th because of a non-ISM absorption feature at 3957. 3 A. 
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Fig. 12. Same as Fig. [8] The spectra around A3951 have ragged baselines, presumably due to instrumental fringing. Only the 
coadded spectrum at A3951 is used for subsequent averaging of R with that from AA232. 



